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A combination of Euler parameter kinematics and Hamiltonian mechanics provides a rigid 
body dynamics model well suited for use in strongly nonlinear problems involving arbitrarily 
large rotations. The model is unconstrained, free of singularities, includes a general poten- 
tial energy function and a minimum set of momentum variables, and takes an explicit state 
space form convenient for numerical implementation. The general formulation may be spe- 
cialized to address particular applications, as illustrated in several three dimensional example 
problems. 


INTRODUCTION 

A variety of rigid body dynamics modeling problems demand consideration of very large 
rotations. Some of the best known examples involve aircraft [1] and spacecraft [2,16], al- 
though the analysis of large rotation dynamics is of generic interest in a wide range of 
applications, including mechanism and machine theory [3,4] and molecular dynamics [12]. 
Most models of rigid body dynamics problems employ Euler angles [5]. Such formulations 
lead to equations of motion which are unconstrained, but which contain singularities [1]. 
Other singular three parameter methods have been developed, including for example those 
of Laning-Bortz-Stuelpnagel [11] and Rodriguez [13]. The presence of singularities in all these 
methods has motivated the development of alternative four parameter modeling schemes [11], 
including Euler parameters [15]. Such formulations replace the three Euler angles with four 

3 Graduate research assistant 

4 Professor 


43 


parameters and an algebraic constraint. This avoids the Euler angle singularities but leads 
nominally to a system level model in differential-algebraic form. 

In an attempt to avoid both singular equations of motion and differential-algebraic sys- 
tems, several authors have presented reformulations of Euler parameter based models, for 
use in three dimensional rigid body dynamics problems. Chang et al. [2], Nikravesh and 
co-workers [7,8,9,10], and Vadali [17] present alternative formulations based on Lagrange’s 
equations. Since Lagrange’s method defines the solution as a path in configuration space, 
and since the Euler parameters are taken as generalized coordinates, this approach starts 
with a differential system of order eight (four second order equations for the rotational dy- 
namics of a single rigid body) augmented with a single algebraic constraint. Nikravesh and 
co-workers begin from this starting point and proceed to find a closed form solution for the 
Lagrange multiplier associated with the algebraic constraint, resulting in an unconstrained 
formulation of order eight. They do not include a potential energy function in the system 
Lagrangian. Similar results are obtained by Vadali. Proceeding in a different manner, Chang 
et al. introduce as quasi-velocity variables the rigid body angular velocities in the body fixed 
frame, and project the original order eight Lagrange equations onto an order seven subspace. 
In the process they eliminate the unknown Lagrange multiplier. 

As an alternative to Lagrange’s equations, a Hamiltonian formulation of rigid body dy- 
namics with Euler parameters has been proposed by Morton [6]. However his final formula- 
tion is of order eight, and includes a superfluous momentum variable as well as a ‘generally 
arbitrary’ unspecified scalar parameter. It appears that no previous work has attempted to 
revise or improve upon the Morton formulation. 

The usefulness of formulations based on Hamilton’s canonical equations is well recognized 
[3]. They offer an explicit state space description of system dynamics problems which is: 
(a) convenient for numerical integration, (b) well suited for coupling to automatic control 
system models, and (c) energy based and hence providing clear physical insight. Recognizing 
these strengths, a revision and extension of existing Hamiltonian formulations for rigid body 
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dynamics is of generic interest. The present paper presents such work, deriving unconstrained 
Hamilton’s equations for the three dimensional dynamics of a rigid body in terms of Euler 
parameters, and hence suitable for use in simulations involving arbitrary rotational motion. 
The derivation avoids any requirement to determine the Lagrange multiplier associated with 
the Euler parameter constraint. No arbitrary parameters are introduced, and the final 
rotational formulation is of order seven. A general potential energy function and nonpotential 
virtual work effects are included in the model. Validation and application of the method is 
illustrated here in several three dimensional example problems. 

KINEMATICS 

This section defines the kinematic variables of interest, and recalls a number of well known 
kinematic relations [1], for use in succeeding sections. 

The position and orientation of an arbitrary rigid body is described here in terms of seven 
generalized coordinates, namely the Cartesian components of the center of mass vector ( c ) 
and a four component vector of Euler parameters ( e ) 

c = [ ci c 2 c 3 ] T , e = [ e 0 ei e 2 e 3 } T (1) 

Knowledge of the Euler parameters determines a (nonunique) set of Euler angles (0, 9, ij>) 
for the body, associated with a 3-1-3 rotation sequence, via the relations 


0 

= tan *( — ) + tan x ( 


(2) 


= fa77, _1 ( — ) — tan~ l { 

A 

(3) 


e 0 

ei 


9 

= 2 sin~ 1 ( \/ ei 2 + e 2 2 

) 

(4) 


The Euler parameters define an orthogonal rotation matrix (R) which relates the vector of 
components (a) of a first order tensor described in a fixed Cartesian coordinate system to a 
corresponding vector of components (a) described in a body fixed co-rotating frame, using 

a = R a (5) 
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where 


with 
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( 6 ) 
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The four Euler parameters are not independent, and satisfy the constraint equation 

e T e = 1 (9) 

which then implies 

G G r = I (10) 

where I is an order three identity matrix. In addition, G and e and their time derivatives 
satisfy the identities 

G e = 0, G e = — G e (11) 

The kinematic equations [1] which relate the time derivatives of the Euler parameters 
to the components of the angular velocity vector ( uj ) of the rigid body, expressed in the 
body-fixed co-rotating frame, are 

w = 2Ge, e = ^ ^ T ^ (12) 

Finally note that the skew-symmetric matrix 12, with axial vector ay which satisfies 

f2v = wxv (13) 

for any vector v, is related to the Euler parameters and their time derivatives by 

n = 2 G G t = — 2GG t (14) 
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The next section defines kinetic and potential energy functions and hence the Hamiltonian 


for the system of interest. 


KINETIC AND POTENTIAL ENERGY 

The complementary kinetic energy for the rigid body may be expressed as 

T* = - m c T c + - uj 1 J lo 
2 2 


(15) 


where m is the mass and J is a constant matrix of components for the moment of inertia 
tensor of the body, referred to the co-rotating frame. In terms of the Euler parameters and 
their time derivatives, 

r = ^mc T c + 2e T G T JGe (16) 

which has the form T* = T* ( c, e, e ). It follows that the generalized momenta are 


dT* 

P = ~ rp = rn c , 

a c 


dT* T 

g= — =4G r JGe 
ae 


(17) 


Note that the identities (10) through (12) require 


g = 2 G t J uj , 


uj = - J^Gg 
2 B 


(18) 


Since the complimentary rotational kinetic energy may also be expressed, using equation 
(11), as 

r; ot = 2e r G r JGe (19) 

then the Euler parameter dependence of T* defines the partial derivative 


dT* 

k = =4G T JGe 

de 

The kinetic energy of the body is defined via the Legendre transform 


( 20 ) 


T = p J c + g T e - T 


so that the preceding results lead to the canonical form T = T ( p, g, e ) which is 

T- ( p'p I ‘ g''G r .l r G g 


( 21 ) 


( 22 ) 
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and require that (see the appendix) 


, dT 


( 23 ) 


For the mechanical systems considered here, the potential energy function has the general 
form 

V — V(c, e) (24) 

and the system Hamiltonian is 

H = T + V (25) 

The next section introduces a virtual work expression, to account for nonpotential effects. 


NONPOTENTIAL VIRTUAL WORK 

The quasi-coordinates q associated with the co-rotating components of the angular veloc- 
ity vector are defined by 

q = u) (26) 

In terms of the latter coordinates, and the center of mass coordinates, the nonpotential 

virtual work due to the imposed forces f (t) and torques T(i) is 

5W nc = f(t) T hc + T(f) T hq (27) 

Since 

hq = 2 G he (28) 

the virtual work expression which defines the generalized nonpotential forces is 

8W nc = f (t) T hc + 2 [G t T (t) ] T he (29) 

Note that damping effects may contribute additional terms to the nonpotential virtual work, 
in which case the nonpotential forces may depend on the generalized coordinates and ve- 
locities. In addition the presence of nonholonomic constraints may introduce terms which 
depend on unknown Langrange multipliers. The last problem discussed in the examples 
section illustrates the effects of both damping and nonholonomic constraints. 


48 


The next section derives Hamilton’s equations for the system. 


HAMILTON’S EQUATIONS 

The system Hamiltonian has the form H = H ( p, c, g, e ) and the canonical Hamilton’s 
equations are 

dH dH 

P = -^ + f p , c = — (30) 


and 


g = 


dc 

dH 


f 9 . 


dp 

dH 


e = 


(31) 


de ’ d g 

where f p and f 9 are the nonpotential generalized forces associated with the virtual work and 
any applied constraints. The Euler parameter constraint has the rate form 


e 1 e = 0 


(32) 


Introducing a Lagrange multiplier A, the latter constraint combines with the virtual work 
expression to yield 


f p = f(t) 

f 9 = 2 G T T(f) + Ae 

so that for the derived Hamiltonian the momentum balance equations are 

dV 

p = -^ + f(<) 

dV 

g = + k + 2G T T(t) + Ae 

ae 


(33) 

(34) 


(35) 

(36) 


The last equation includes an unknown Lagrange multiplier and a superfluous momentum 
variable. These variables are eliminated by introducing the three-momentum vector 


h = Jw 


(37) 


whose time derivative is 


h= ^Gg + ^Gg 


(38) 
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With equations (11) and (36) this yields 



( 39 ) 


or 


11 1 dV 

h = --n h + - Gk + T (f)--c;^ 


(40) 


which eliminates both A and the four component momentum vector. 


The final unconstrained Hamiltonian model is 


P 


h 



(41) 


(42) 


c 



(43) 


e 




(44) 


Given a potential function and the virtual work, the preceding explicit equations may be 
integrated to simulate the system response. 

As outlined in the introduction, the rigid body dynamics formulation derived here com- 


nian and Lagrangian methods generally simplify the model formulation process, in particular 
when geometric nonlinearities are important. Such nonlinearities arise for example when 


normally most convenient for numerical integration. Euler parameters offer a singularity 
free description of rotational displacements, and are therefore preferred over Euler angle 
models in large rotation applications. The cost is of course the need to integrate an addi- 
tional state equation for each rigid body, since the Euler parameters are a quaternion. The 
present combination of Hamiltonian mechanics and Euler parameter kinematics is therefore 
of most interest in the formulation and numerical integration of models for strongly nonlin- 
ear mechanical systems. The model developed here is unique in its combination of features: 


bines the advantages of Hamiltonian mechanics and Euler parameter kinematics. Hamilto- 


large rotations or hyperelastic devices are of interest. As compared to Lagrangian meth- 
ods, Hamiltonian methods offer an explicit state space description of the system dynamics, 
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unconstrained, free of singularities, incorporating a general potential energy function, em- 
ploying a minimum set of momentum variables, and taking an explicit state space form 
convenient for numerical implementation. 

EXAMPLE PROBLEMS 

Application of the Hamiltonian formulation developed here is illustrated in four examples. 
The first example compares an Euler angle based model of a rotating disk problem to the 
present Hamiltonian formulation, both to validate the present approach and to illustrate 
in a simple case an Euler parameter description of rotational displacement. The second 
example solves a classical rigid body dynamics problem for the three dimensional motion of 
a torque-free body, for comparison to the published numerical solution of Morton [6] and to 
the partial analytical solution of Thompson [16]. This problem again validates the present 
approach, and compares a discontinuous Euler angle description of three dimensional rigid 
body motion to the continuous Euler parameter characterization adopted here. The third 
example models a spinning top in a uniform gravitational field, a problem described in many 
advanced dynamics tests, and validates the present formulation in a three dimensional rigid 
body motion involving both translation and rotation. Here a numerical solution for the last 
cited problem is obtained using a time step identical to that employed by Simo and Wong 
[14], but without resort to their symplectic integration algorithm. The last example considers 
a problem of practical importance in the design of gyroscopic seekers, namely the motion of 
a freely precessing body with a viscous ring nutation damper [2], This example calls for the 
application of nonholonomic constraints. Here we develop an explicit state space model, as 
compared to the implicit Lagrangian formulation [2] of Chang et al. 

The first example models the free vibration of a rigid circular disk of radius r, rotating 
about a fixed point, and attached to a linear spring of stiffness k (see Figure 1). The potential 
energy function is 

V = ^ ky 2 p = l,kr 2 [ 2eie 2 + 2e 0 e 3 ] 2 (45) 
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where y v denotes the vertical displacement of the point of attachment of the spring, measured 
in a Cartesian coordinate system whose origin lies at the center of the disk. Note that the 
indicated Euler parameter dependence of the potential energy is obtained using equation 
(5). Assuming the model parameters and initial conditions listed in Table 1, the motion was 
simulated by integration of a Newtonian model based on the Euler angle cj) 

J(j)+^kr sin{ 20) = 0 (46) 

and by integration of the Hamiltonian relations (42) and (44). The two computed results for 
the angular momentum are compared in Figure 2, showing excellent agreement of the Euler 
angle based and Euler parameter based solutions. Figure 3 shows the computed variation of 
the Euler parameters with time. 

The second example models the torque free motion of a rigid body, for the inertial prop- 
erties and initial conditions listed in Table 2. A partial analytical solution for this classic 
problem is known and can be expressed in terms of elliptic functions [6,16]. Figures 4 and 
5 shown the time variation of the angular momenta and Euler parameters computed using 
the present Hamiltonian formulation. Figure 6 plots the implied Euler angles, emphasizing 
the discontinuous nature of the latter variables. Table 3 shows excellent agreement of the 
analytical and numerical solutions for the amplitudes and periods of the angular momenta, 
in three dimensional motion. 

The third example models the translational and rotational motion of a spinning top in 
a uniform gravitational held. This problem is described in many advanced dynamics texts 
[4], and is used by Simo and Wong [14] to evaluate their symplectic numerical integration 
scheme. The potential energy for the system is 

V = W z c (47) 

where W is the weight of the top and z c is the vertical coordinate of the center of mass. The 
simulation parameters and initial conditions for the problem are listed in Table 4. Figures 7 
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and 8 show numerical results for the normalized angular momentum components and center 
of mass coordinates, obtained by integration of Hamilton’s equations (41) through (44), 
using a fourth order Runge-Kutta method. Table 5 compares an approximate analytical 
estimate of the nutation and precession frequencies for the top, provided by Goldstein [4], 
to the present numerical results. The present numerical results are identical to those plotted 
by Simo and Wong [14], and are obtained using the same time step, but without resort to 
their symplectic integration scheme. 

The fourth example considers the rotational motion of a rigid rotor damped by a partially 
filled mercury ring damper. The reader is referred to Chang et al. [2] for a detailed discussion 
of this problem, and its application in the analysis of gyroscopic seekers. We focus here on 
the formulation of a dynamic model for the system analyzed in reference [2], The paragraphs 
which follow develop an explicit Hamiltonian model for this system, an alternative to the 
implicit Lagrangian model of Chang et ah, adopting their stipulated assumptions on stored 
energy functions, energy dissipation, and kinematic constraints. 

The rotor is modeled as a rigid circular cylinder with a fixed center of mass located at 
the origin of a global XYZ coordinate system. The partially filled mercury ring damper 
is a cylinder of mean radius R , co-axial with and external to the rotor, with a centroid 
displaced a distance L along the rotor axis from the rotor center of mass location. Body- 
fixed coordinate systems for the rotor (xyz) and ring (uvz) are co-located at the centroid of 
the damper, where the z direction is aligned with the rotor axis. The partial mercury ring is 
free to rotate about the rotor axis, subject to a damping torque which is linear in the axial 
angular velocity difference between the rotor and the ring, but is otherwise constrained to 
move with the rotor. Hence the orientations of the body-fixed axes systems which co-rotate 
with the rotor and the ring differ only by an angle (3, which describes the axial rotation of 
the ring with respect to the rotor. 
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The assumed complimentary kinetic energy for the system is [2] 


r = \ ^ T J T « + \ (48) 

where U3 and are angular velocities for the rotor and ring and 
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are constant moment of inertia matrices for the rotor and ring. All four quantities are 
described in the respective rotor and ring body-fixed co-rotating coordinate systems. 

The assumed potential energy for the system, due to the gravitational potential of the 
mercury, is [2] 

V = -m g c R T B T r c , r c = [R sin{^/2)/{^/2), 0, L} T (50) 


where m is the mass of the mercury, g c is a constant gravity acceleration vector described 
in the fixed XYZ frame, R is the rotation matrix of equation (6), whose Euler parameters 
(e) refer to the rotor-fixed frame, r c is a constant vector which locates the mercury center 
of mass, 7 is the angle which subtends the mercury arc (symmetric about the u axis), and 
B is an orthogonal matrix which defines the transformation of vector components from the 
rotor-fixed to the ring-fixed frame 


B 


cos(/3 ) sin(fl) 0 
— sin(/3 ) cos(fi) 0 

0 0 1 


(51) 


Note that V = V(e,/3). The virtual work for the system, due to damping at the ring-rotor 
interface, is [2] 

5W = -C d R 2 ft 6P (52) 


where Cd is an empirical dimensionless damping coefficient. 

Given the preceding modeling assumptions, Hamilton’s equations for the rotor and ring 
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system are 


h = 

-fih-^G^ + T 

2 oe 

(53) 

h m 

+ T m 

(54) 

e = 


(55) 

0 = 

dV m 
~~dp + p 

(56) 


where h and h m are angular momenta for the rotor and ring 

h — J uj , h m = 3 m uj m (5 7) 

and T. T m , and Tp are nonpotential forces due to damping and kinematic constraints. Note 
that the degenerate form of Hamilton’s equation for (3 is due to the fact that the latter 
generalized coordinate, which appears in the potential energy function, is not associated 
with a corresponding generalized momentum variable. The kinematic constraints are [2] 

Umi = cos{(3) uii + sin(P) lu 2 , Cc ; m2 = — sin((3 ) LUi + cos{(3 ) u 2 , (58) 

and 

P = U>m3 - ^3 (59) 

They quantify the aforementioned modeling assumption that the ring moves relative to the 
rotor only in axial rotation. 

An explicit state space model may be obtained by application of the constraints, as follows. 
Introducing a Lagrange multiplier /i for the constraint (59), and accounting for the virtual 
work, requires 

T /3 = n - C d R 2 (3, T = n c, T m = -fx c (60) 

where c denotes the vector [0,0, 1] T . The degenerate Hamilton’s equation for (3 therefore 
determines the Lagrange multiplier as 

8V 

H= C d R 2 (3+ — (61) 


55 


Hamilton’s equation (54) for the ring angular momentum may now be written in the form 


<jJ r 


— J™ ^m.Jr 




h Jr 


(62) 


Since the constraints specify both (3 and the first two components of the ring angular velocity 
vector as functions of the set (f3,coi,U2,^m 3 ), the third of equations (62) is an evolution 
relation for the unknown uJm 3 . Combining the third of equations (62) with the constraint 
equation (59), the constitutive relations (57), and Hamilton’s equations (53) and (55), the 
result is an explicit state space model of order nine for the ring-rotor system. The final state 
equations are 


h = 


e = 


P = 


^m3 — 


_f!h_ 2 G ^ + 


C, U 1 


h 3 \ dV 

“ ” 3 ~TJ + W. 


-C T J- lh 

h'i 

W m 3 V~ 

"3 

^ (j m 


C d R 2 


h 3 \ dV' 


c T Jjc 


where 


f2 — fi(h), G — G(e), Q, m — Q m (u: m ) 


(63) 

(64) 

(65) 

( 66 ) 

(67) 


with 


cVmi = cos({3) ^ + sin(P) 

J 1 J 2 


w m2 = - sin(/3 ) ^ + cos(f3) (68) 

J 1 J 2 


Note that implicit model of reference [2] is also of order nine but employs a different set of 
state variables. 


CONCLUSION 

The present paper has derived and applied a new Hamiltonian formulation of rigid body 
dynamics problems, based on Euler parameter kinematics. Euler parameter kinematics pro- 
vide a singularity free description of three dimensional rigid body motion, accommodating 
arbitrarily large rotations. When combined with Hamiltonian mechanics, the result is an 
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energy based modeling approach well suited to address problems with complex geometric 
nonlinearities. As compared to previous work, the formulation derived here offers a unique 
combination of features. It avoids the introduction of algebraic constraints and unspecified 
parameters, includes a general potential energy function, incorporates a minimum set of 
momentum variables, and takes an explicit state space form convenient for use in control 
related applications. 
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APPENDIX 


A complimentary kinetic energy expression (T*) with the functional form 


T* = T*(e,f), f = e 


( 69 ) 


has the total differential 


f)T* T 

dT* = g T di + — de , 
de 


g = 


dT* 

w 


(70) 


The corresponding kinetic energy function (T) is determined by the Legendre transform 


T = g T f - T* 


(71) 


and has a total differential defined by 


BT* t 

dT = g T di + f r dg — dT* = f r dg — de (72) 

oe 


as well as the canonical form 


rrTi 3 T t , 3T t , 
dT = — dg + — de 
a g oe 


It follows that 


dT dT _ dT* 

dg ’ de de 


(73) 


(74) 
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Parameter 

value 

Mass moment of inertia (kg m 2 ) 

J = 2 

Radius of the disk (m) 

r — 1 

Stiffness of the spring ( N/m ) 

k = 10 

Initial displacement ( degrees ) 

0 = 30 

Initial momentum (kg m 2 rad/s) 

h = 0 


Table 1: Model parameters and initial conditions for the first example problem 


Parameter 

value 

Mass moments of inertia (kg m 2 ) 
Initial Euler parameters 
Initial momenta (kg m 2 rad/ s) 

J\ = 400, J 2 = 307.808385, J 3 = 200 
e 0 = 1, ei = e 2 = e 3 = 0 
h i = 346.4101616, h 2 = 0, h 3 = -200 


Table 2: Model parameters and initial conditions for the second example problem 


Variable 

exact solution 

numerical solution 

Magnitude of h\ (kg m 2 rad/s) 

346.4102 

346.38 

Magnitude of h 2 (kg m 2 rad/ s) 

365.447 

365.44 

Magnitude of h 3 (kg m 2 rad/s) 

200.0 

199.975 

Period of hi (s) 

9.3393 

9.35 

Period of h 2 (s) 

18.6786 

18.68 

Period of h 2 (s) 

18.6786 

18.68 

Minimum of h 3 (kg m 2 rad/ s) 

162.6296 

162.6342 


Table 3: Exact versus numerical results for the second example problem 


61 


Parameter 

value 

Weight ( kgm/s 2 ) 

Mass moment of inertia (kg m 2 ) 
Initial Euler parameters 
Initial angular momenta (kg m 2 rad/s) 

W = 20 

Ji = 5, J 2 = 5, J 3 = 1 

e 0 = cos(0.15), e\ = sin(0.15), e 3 = 0, e 4 = 0 
hi = 0, h 2 = 0, h 3 = 50 


Table 4: Simulation parameters and initial conditions for the third example problem 


Variable 

analytical approximation 

numerical simulation 

Nutation frequency (rad/s) 

10.00 

9.24 

Precession frequency (rad/ s) 

0.40 

0.4136 


Table 5: Approximate analytical versus numerical results for the third example problem 
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LIST OF FIGURES 


Figure 1. First example problem, rotating disk with a translational spring suspension 

Figure 2. First example problem, comparison of Euler angle based and Hamiltonian solu- 
tions for the angular momenta versus time; the computed Hamiltonian component h z agrees 
with the Euler angle solution (diamond symbols) obtained using equation (46), while the 
Hamiltonian components h x and h y are identically zero 

Figure 3. First example problem, numerical solution for the Euler parameters versus time; 
note that the Euler parameters e± and e2 are identically zero, and along with the computed 
solutions for eo and e 3 determine the nonzero Euler angle in accordance with equation (2) 

Figure 4. Second example problem, torque free motion of a rigid body, numerical solution 
for the angular momenta versus time 

Figure 5. Second example problem, torque free motion of a rigid body, numerical solution for 
the Euler parameters versus time; note that the Euler parameters are continuous functions 

Figure 6. Second example problem, torque free motion of a rigid body, computed Euler 
angles versus time; note the discontinuities in the Euler angles 

Figure 7. Third example problem, translation and rotation of a spinning top, numerical 
solution for the normalized components of the angular momentum versus time 

Figure 8. Third example problem, translation and rotation of a spinning top, numerical 
solution for the center of mass position versus time 


63 


Figure 1: First example problem, rotating disk with a translational spring suspension 



time in secs 

Figure 2: First example problem, comparison of Euler angle based and Hamiltonian solu- 
tions for the angular momenta versus time; the computed Hamiltonian component h z agrees 
with the Euler angle solution (diamond symbols) obtained using equation (46), while the 
Hamiltonian components h x and h y are identically zero 
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time in secs 


Figure 3: First example problem, numerical solution for the Euler parameters versus time; 
note that the Euler parameters e± and e 2 are identically zero, and along with the computed 
solutions for eo and e 3 determine the nonzero Euler angle in accordance with equation (2) 


66 


400 



Figure 4: Second example problem, torque free motion of a rigid body, numerical solution 
for the angular momenta versus time 
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Figure 5: Second example problem, torque free motion of a rigid body, numerical solution for 
the Euler parameters versus time; note that the Euler parameters are continuous functions 
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Figure 6: Second example problem, torque free motion of a rigid body, computed Euler 
angles versus time; note the discontinuities in the Euler angles 
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Figure 7: Third example problem, translation and rotation of a spinning top, numerical 
solution for the normalized components of the angular momentum versus time 
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Center of mass location 



Figure 8: Third example problem, translation and rotation of a spinning top, numerical 
solution for the center of mass position versus time 
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